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ABSTRACT 


The  concept  of  shock  wsve  physics  in  solids  is  generalized  to 
include  thermal  (second  sound)  shocks,  electromagnetic  shocks,  as  well 
as  the  more  generally  understood  mass  density  shocks.  It  ir  speculated 
that  thermal  shocks  may  be  made  to  occur  in  room  temperature  metals  by 
causing  a  high  power  pulsed  laser  to  impact  on  the  metal  surface.  A  laser 
pulse  width  of  lO"1^  seconds  is  derived  for  the  effect  to  occur. 

The  state-of-the-art  with  respect  to  understanding  the  co>  ipaction 
behavior  of  porous  materials,  the  microscopic  explanation  of  electro¬ 
mechanical  effects  (shock  polarization)  in  a  variety  of  media,  and  the  state 
of  our  understanding  of  the  physics  of  the  detonation  process  is  discussed. 


L  INTRODUCTION 


The  decades  following  the  outbreak  of  the  second  work-  war  have 
seen  the  subject  of  shock  wave  physics  become  more  than  just  u  large 
amplitude  step- child  of  hydrodynamics.  Shock  wave  phenomena  are  now 
known  to  be  important  to  solar  physics1,  astrophysics0,  possible  controlled 
fusion  mechanisms'1* 4  ,  solid-state  physics5*®,  the  physics  of  explosives^* 
supersonic  flight,  and  a  host  of  modern  military  applications1. 

Prior  to  the  second  world  war  shock  wave  effects  vere  treated  in 
a  typical  course10  on  hydrodynamics  "merely  tor  the  rea-  on  that  all  actual 
fluids  are  more  or  less  compressible".  We  can  appreciate  the  growth  of 
the  subject  of  shock  wave  physics  by  realizing  that  whe»  e  many  phenomena 
are  concerned  •-*•<,>  consider  linear  physics  merely  because  it  represents  the 
limiting  ease  of  what  we  really  wish  to  s  tudy.  Further,  *  'here  hydrodynamics 
dealt  exclusively  with  mass  flow,  present  shock  wave  physics  is  much  more 
general;  we  shall  be  concerned  also  with  shock  waves  aosociated  with  the 
propagation  of  temperature  waves*0  (in  the  absence  of  a  mass  disturbance), 
and  shock  waves  associated  with  the  propagation  of  electromagnetic  waves10 
in  ?  medium  of  field  permittivity  and  susceptibility.  In  addition  hydrodynamics 
as  a  subject  does  not  consider  phenomena  on  a  microscopic  level  as  it  is  a 
combination  of  continuum  mechanics  plus  thermodynamics.  In  shock  wave 
physics  we  have  bee  me  very  much  aware  of  the  role  of  microscopic  mechanisms 
on  the  observed  macroscopic  observables!**. 


IL  A  UNIFIED  APPROACH 


That  a  Propagating  disturbance  takes  the  form  of  a  shock  is  a 
consequence  of  the  equation  ot  state  appropriate  to  the  material  in  question 
and  the  type  of  distur-hancV  being  studied.  In  general,  from  a  systems 
point  of  view,  a  nonljk response-stimulii  relation  can  give  rise  to  shock 
like  phenomena.  * 

* 

For  a  mess  «)  asity  change  propagating  into  a  fluid,  the  shock 
velocity  Us  is  given 


.A  /C.  /  p  -  p.  \ 

V  U  -/  li  (  f  1 

s-/p f  Vvpi/' 


while  the  propagation  velocity  for  an  infinitesmal  disturbance  with  respect  to 
any  local  equilibrium  state  of  the  material  is  given  by  the  "sound  speed" 


sound  speed  = 


In  equations  (1)  and  (2),  p  is  density,  P  is  pressure,  and  the  subscripts  i  and  f 
refer  to  before  and  after  the  shock,  respectively.  If  in  the  above  we  think  of 
AP  as  being  the  stimulus  and  A p  as  the  response, then,  before  further  detailed 
consideration  of  shock  propagation,  it  seems  worthwhile  to  generalize  equations 

(1)  and  (2)  to  other  systems.  Table  I  demonstrates  the  generalization. 

TABLE  I 

GENERALIZATION  OF  STIMULUS  AND  RESPONSE 


Stimulus 


Response 


(1)  Pressure  Change,  AP  Density  Change,  Ap 

(2)  Temperature  Change,  AT  Entropy  Change,  AS 

(3)  Electric  and  Magnetic  Changes  in  the  induced  fields 

Field  Changes,  AE  and  AH  A’O  and  A  R 


i 


Where  the  three  systems  given  in  Table  I  are  these  which  will  be  considered 
in  some  detail  here,  they  do  not  exhaust  the  relevant  possibilities;  a  non¬ 
linear  current-voltage  relationship  can  give  rise  to  shock  effects  which  are 
of  practical  importance  to  transmission,  line  theory1®  or  D.  C.  to  A.  C. 
converters  utilizing  semiconductor  materials*®*  17. 


Based  upon  equations  (1)  and  (2),  and  Table  I,  we  might  expect  the 
local  "sound"  velocities  and  shock  velocities  to  be  given  by: 


Souadvclocity18  “y 

/  b  T  ' 

.  «  t 

(3) 

d  S 

Shock  velocity  o<  j 

r  t7  -  t"-1 
f  i 

(4) 

/  f  Jt  . 

J  sf  -  S{ 


Electromagnetic  Propagation: 


Sound  velocity19  , 

J  bD  bB 


Shock  velocity 


(5) 


In  equations  (3)  to  (8)  inclusive,  «  means  proportional  to;  the  generalization 
leading  to  the  equations  does  not  guaranty  that  numerical  constants  or  other 
factors  will  not  appear  in  the  exact  equations* 

To  the  large  majority  of  those  involved  in  shock  wave  physics,  the 
inclusion  of  phenomena  other  than  that  associated  with  a  propagating  mass 
discontinuity  requires  some  explanation.  There  are  three  strong  reasons  tor 
the  generalization  we  are  making,  and  they  arc  given  below: 


a.  Both  the  thermal  and  electromagnetic  shocks  get  coupled 
to  the  mass  density,  so  that  the  thermal  and  electromagnetic  shocks  can 
give  rise  to  the  more  famila .  mass  density  discontinuity.  The  thermal 
expansion  coefficient  is  a  coupling  coefficient  for  the  thermal  shock,  while 
charge  density  (which  is  made  up  of  ions  and  electrons),  current  density, 
and  bulk  polarization  couple  she  field  terms  in  Maxwell's  equations  with  the 
mass  flow. 


b.  A  modern  method  of  introducing  a  propagating  mass  dis¬ 
continuity  is  io  irradiate  a  target  material  with  a  pulse  of  high  energy 
electrons2®  ,  or  a  pulse  of  photons  (e.  g. ,  from  a  laser  beam).  We  will 
predict  in  this  report  thai  a  sharp  discontinuity  in  thermal  energy  deposition, 
for  example  with  skin  depth  rise  time  as  occurs  with  a  laser  pulse,  can  give 
rise  to  a  temperature  (entropy)  shock  in  addition  to  usual  mass  discontinuity 
which  results  from  thermal  pressure.  Pressure  is  equivalent  to  energy 
density.  Since  the  mass  discontinuity  and  thermal  discontinuity  propagate 

at  different  velocities,  the  neglect  of  the  temperature  shock  could  cause  an 
error  in  estimating  transport  associated  with  the  mass  discontinuity. 

c.  The  formalism  for  the  three  different  shock  phenomena 
listed  in  Table  I  are  essentially  the  same,  even  to  employing  characteristics 
in  seeking  mathematical  solutions13,  21.  In  addition,  the  phenomena  are 
essentially  the  same  although  particular  mechanisms  are  different;  jump 
conditions,  viscosity,  and  nonlinear  equations  of  state  are  all  applicable  to 
each  phenomena. 


The  Electromagnetic  Shock  -  Let  us  consider  the  c,wo  time  varying  Maxwell's 

equations  for  electromagnetic  field  propagation 
in  a  medium. 


V  x  H  = 


a  d 

d  t 


(7) 


V  X  E 


1  <LB 
c  ,3  t 


(8) 


22 

where  gaussian  units  have  been  used*'  ,  e  is  the  velocity  of  light  in  vacuum, 

— ^  — x  ■  ^ 

J  is  current  density,  E  is  the  electric  field  vector,  H  the  magnetic  field  vector, 
D  and  B  are  the  corresponding  ind*  ced  vectoi  fields,  V  is  the  gradient  operator, 

■  ^  -A 

and  (7  x  H)  denotes  the  usual  "curl  of  H". 


We  now  restrict  ourselves  to  the  one  dimensional  strain  type  of 
problem,  add  allow  for  spatial  variation  only  in  the  Z  (or  Xg)  direction* 

Further  restricting  ourselves  to  the  special  case  J  =  0  gives 


9H2 

i 

3D1 

— and 

aHi_  i 

aD, 

c 

dt 

at 

5^2 

1 

fiBi 

— i.  and 

a  Ej  i 

-  3  —  - 

aB, 

ax3 

c 

at 

ax3  c 

at 

'9) 


(10) 


Here  the  subscripts  1,  2,  and  3  refer  to  components  along  the  X,  Y,  and  Z 
directions,  respectively-  If  we  assume  a  steady  propagating  shock,  with  the 
velocity  Ue  ,  then  we  can  write  dx3  =  Ug  ui  .  With  respect  to  an  observor  on 
the  propagating  front  the  first  equation  of  equation  (9)  and  the  second  equation 
of  equation  (10)  become,  respectively, 


dH2 

u„ 

s--  dD,  , 

C  1 

(11) 

dV 

~  dB2  • 

(12) 

Integrating  equations  (11)  and  (12)  acrons  the  shock  front  gives 


• 

:  2a. 

c 

\°A  • 

(13) 

el 

1 

«jwC: 

iB2)  * 

xl4) 

where  {  j  indicates  the  jump  (i,  «.  ,  change)  in  the  qua  tty  E,  across 
the  3bock  front  Combining  equations  (13)  and  (14)  yields 


(15) 


For  the  propagation  velocity  of  small  amplitude  disturbances  in  Ej  and  E^ 

within  the  larg .)  field  region  where  e  and  fj,  are  to  be  considered  constants, 
we  can  form  the  wave  equation 


1  _ 


a  x. 


s2e] 

at  2 


(16) 


Equation  (16)  follows  directly  from  operating  on  equations  (9)  and  (10).  Thus, 
sound  speed  =  =  c  /d-^1  dHo~  1  _ 

J  3D1  SB2 


(17) 


We  thus  see  that  the  formalism  for  the  electromagnetic  shock  corresponds, 
as  expected,  to  the  formalism  for  a  propagating  mass  discontinuity  as 
expressed  in  equations  (1)  and  (2). 

If  we  now  include  a  nonvanishing  current  density,  then  equation  (11) 

becomes 


-dH  = 

2  c 


4tj 


J]  dx3 


U, 


dD 


1  > 


(18) 


while  equation  (12)  remains  unchanged.  The  shock  velocity  thus  becomes 


(19) 


where  the  integration  is  carried  out  over  the  shock  front  to  the  points  where 
the  jump  values  are  observed.  It  is  not  at  all  surprising  that  enters  into 

equation  (19)  as  it  does;  currents  are  always  associated  with  their  surrounding 
magnetic  fields.  In  evaluating  equation  (19)  particular  model.*:  must  be  con¬ 
sidered  as  the  current  density  is  dependent  upon  the  field  variables. 


The  current  density  plays  the  same  role  in  the  electromagnetic 
shock  that  the  flow  gradient  viscosity  plays  in  a  mass  shock,  namely  the 
sn.oothiag  of  the  discontinuity.  The  current  density,  through  the  electrical 
conductivity  of  the  material,  gives  rise  to  a  viscosity  and  associated 
penetration  depth2®.  Thus,  the  presence  of  current  density  competes  against 
the  nonlinearity  in  e  and  p  in  establishing  the  thickness  of  the  shock  fiont 
Thus,  the  picture  is  a  large  amplitude  plane  wave  electromagnetic  pulse 
(for  example  a  step  function)  is  by  some  means  caused  to  impinge  upon  a 
solid  and  propagate  through  it.  Inertial  effects  are  coup'ed  to  the  electro- 
magnetic  pulse  via  the  ion  and  electron  current  densities  which  make  up  Jj  . 
After  the  pulse  has  propagated  into  the  medium  of  finite  distance  we  expect 
a  steady  stats  to  be  reached  with  the  rise  time  of  the  shock  being  controlled 
as  mentioned  ahove. 


The  Second  Sound  Shock  -  Second  sound  is  typically  a  low  temperature 
phenomena,  and  is  usually  associated  with  the  suporfluid  stat.o  of  liquid  helium. 
In  the  two  fluid  models  for  liquid  helium  the  mass  density  is  broken  down  mto 
"normal"  and  "superfluid "  parts 


P  =  Ps 


(20) 


where  the  subscripts  s  and  n  refer  to  the  two  parts  mentioned  above.  It  is 
obvious  from  equation  (20)  thatpg  and  p^  can  be  simultaneously  varied  in  such 

a  way  that  keeps  the  total  mass  density  constant  The  physical  mechanism 

for  vaiying  the  ratio  p  is  temperature  change,  with  p  — *■  0  as  T  — *■  0  °K. 

8  D 


Figure  1  gives  the  temperature  dependence  of  p 


8 


and  p  . 
n 


i 
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Figure  1  -  Tbe  temperature 
dependence  of  and 

pn  for  liquid  helium. 
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In  the  two  fluid  models  for  liquid  helium  the  superfluid  component 
has  ^ero  entropy,  while  the  normal  component  has  the  usual  thermal 
excitations  of  a  dense  gas  or  liquid  with  accompanying  disorder  and  entropy. 

Since  (p  +  pQ  )/p  in  figure  1  has  the  property  of  equaling  unity  (i.  e. ,  one) 
a  temperature  change  stimulus  will  cause  an  entropy  change  response  without 
any  corresponding  change  in  mass  density. 

Under  the  more  everyday  conditions  of  heat  flow  a  temperature 
gradient  causes  no  propagating  entropy  response,  and  heat  flow  takes  place 
instantaneously  (L  e. ,  with  an  infinite  propagation  velocity?)  as  described 
by  the  Fourier  equation  of  heat  conduction,  ^or  liquid  helium  below  the 
critical  temperature  of  2. 17° K,  however,  a  propagating  entropy  response 
exists.  Such  a  phenomenon  was  first  observed  by  Peshkov24.  The  difference 
between  the  liquid  helium  case  and  the  more  everyday  heat  conduction  by 
diffusion  case  is  found  in  the  difference  between  the  thermal  diffusivities  (this 
point  is  expanded  upon  later  in  this  report).  In  turn,  the  difference  in 
diffusivities  is  related  to  the  difference  in  specific  heats  and  thermal  con¬ 
ductivities. 

Considerable  experimental  and  theoretical  work  has  been  carried  out 
since  the  original  Peshkov  experiments,  and  it  is  now  an  everyday  occurrence 
for  a  second  sound  paper  to  appear  in  Physical  Review  Letters,  the  Physical 
Review,  or  a  similar  journal.  The  detailed  physics12,  25  of  second  sound, 
second  sound  attenuation,  and  second  sound  shock  has  long  since  been  considered. 

I’d  lowest  order  the  conservation  of  momentum  combined  with  the 
conservation  of  entropy  can  be  shown12,  2^  to  lead  to  the  wave  like  equation 


where 


(21b) 


2 

For  an  ideal  second  sound  Iv  p  =  0)  disturbance  propagating  in  the  z  direction, 
and  without  any  variation  in  the  x  and  y  directions,  equation  (21a)  becomes 


i&.  =Ps  g2  /  9T  \  d2S  . 

at2  Vas /f  az2 


(22) 


One  thus  predicts  a  second  sound  velocity,  Ug.of 


(23) 


in  agreement  with  the  general  ideas  which  weDt  into  equation  (3).  Early 
theoretical  predictions  by  Landau  (1941)  give?*  12>  25  a  ]ow  temperature  limit  of  U 
of 


limit  U  =  Uj 
T— *0°K  2  JV~  } 


(24) 


where  is  the  valocity  of  ordinary  or  first  sound  (i.  e. ,  the  velocity  of 
equation  (2)  at  constant  entropy).  It  is  now  known*2,  25  that  for  Ts  0.  5°K,  U2 
actually  rises  above  the  value  given  by  equation  (u4)  and  may  be  due  to  moan 
free  path  effects. 

To  lowest  order  when  the  conservation  of  entropy  and  the  conservation 
of  momentum  are  combined  with  a  pure  phonon2®  field  (i.  e, ,  no  anharmonic 
effects  or  other  excitations  such  as  rotons2^)  one  arrives12  at  a  sound  velocity 
given  by  equation  (24).  Anharmonic  effects  are,  however,  expected  to  be 
present  in  liquid  helium  at  any  finite  temperature  (after  all  ’t  is  a  liquid  the 
nucleii  of  which  are  not  iu  definite  equilibrium  positions).  Uoton  states  are 
also  occupied.  The  anharmonic  effects  plus  excitations  such  a.,  rotons  leads 
to  deviations  in  the  relevant  equation  of  state  (i.  e. ,  "  vs  S).  These  deviations 
take  the  form  of  nonlinearities  which  in  turn  lead  to  an  amplitude  (i.  e. , 
temperature  change)  dependent  second  sound  velocity  through  equation  (23). 


W'i 
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The  effect  of  the  above  mentioned  nonlinearity  on  a  propagating 
thermal  pulse  for  temperatures  near  but  below  and  above  1.  9°K  is  shown 
in  figures  2a  and  2b,  respectively. 


2a)  Below  1.  9°K: 

LEADING  EDGE 

2b)Above  1.  9°  K: 


Figure  2:  The  propagation  of  second  sound  pulses  of  large  amplitude 
in  liquid  helium  (a)  below  1.  9°K,  and  (b)  above  1.  9°K. 

From  Atkins12. 

The  shock  effect  shown  in  figure  2  is  based  upon  experiments  by  Osborne2®. 

The  situation  depicted  for  below  1<  9°  K  would  seem  to  correspond  to  the  usual 
case  of  a  propagating  mass  discontinuity  where  a  compressive  wave  front 
shocks  up.  The  above  i.  9°K  situation  would  seem  to  correspond  to  the 
infrequently  observed  propagating  mass  discontinuity  case  where  a  rare¬ 
faction  wave  front  shocks  up.  The  cases  (a)  and  fo)  of  figure  2  require, 
respectively 

>0  and  ^5-  <0, 

ds2  a  s2 

which  correspond  to  a_V_  >  o  and  9" V  <  0  t  respectively  for  compressive  and 

Sp2  dp2 

rarefaction  mass  discontinuity  shocks2®. 

We  have  seen  that  second  sound  and  second  sound  shocks  exist  in 
liquid  helium.  Our  more  general  arguments,  which  are  illustrated  in  Table  I, 
would  indicate  that  more  everyday  materials  might  as  well  be  expected  to 
exhibit  the  phenomena.  Chester^  predicts  second  sound  in  the  more  everyday 
materials  providing  that  relaxation  time  and  applied  frequency  conditions  are 
fulfilled.  Chester  notes  that  the  Fourier  heat  equation  is  really  an  approximation 
to 

r  tq  -  K^T  ,  (25) 

at 
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where  q  is  the  thermal  current  density,  ris  a  relaxation  time,  aid  K  the 
coefficient  for  thermal  conductivity.  When  equation  (25)  is  combined  with 
the  continuity  equation  for  heat  transport,  below, 

q  5T  +  ^  •  q  =  0 
St 

then  we  arrive  at  the  modified  wave  equation 

q2T  ,  1  ST  .  ..  1L  v2T  -  0 
at2  t  aT  at 


(26) 


(27) 


In  the  above, C  is  the  heat  capacity  per  unit  volume.  I  .ora  equation  (25)  we 
see  that  in  the  limit  r-»0  we  arrive  at  the  usual  Fourier  heat  equation. 


We  would  now  like  to  explain  the  physics  of  equation  (25).  Consider 
the  temperature  distribution  shown  in  figure  3  below.  Plane  geometry  is 
assumed. 


Figure  3 


fluxes.  Center  region  is  the  region  of  observation. 


Allow  a  heat  flux  meter  tu  be  set  up  somewhere  between  the  temperature  planes 
Tg  and  Tg  .  .  urther,  if  j2  is  the  transporting  particles  (or  quasi-  particles) 

mean  free  path,  let  (Ax)  ~  jj  •  We  ask  ourselves  what  the  flux  meter  measures 
as  a  function  of  time  if  a  switch  is  thrown  at  time  t  =  0  which  changes  the 
temperature  distribution  from  a  homogeneous  one  to  in  winch 


T!  "  T2 
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and  further  that  the  new  temperature  distribution  is  maintained  for  all 
time.  At  time  t  =  0  the  meter  measures  q2g  which  is  the  result  of  the 
local  temperature  distribution  (T2  -  Tg)/(Ax).  If  c  is  the  velocity  of  the 
transporting  particles  then  the  flux  meter  measures,  in  addition  to  q22  , 
flux  contributions  from  regions  I  and  HI  after  a  t4me  t~  Vc.  For  times 
t  >  £/c  the  flux  must  approach  a  steady  state  constant  value  for  a  constan* 
temperature  distribution.  Clearly  then  a  nonconstent  temperature  gradient 
gives  rise  to  a  time  varying  flux.  We  thus  see  that  t  =  0  in  equati  n  (25) 
will  not  explain  the  above  expected  observations. 

In  the  ab'we  paragraph  we  have  seen  tha,  the  flux  contribution,  Aq, 
from  regions  I  and  III  is  of  the  form 

Aq  <*  i  ^ 

dx^ 

For  transporting  velocity  c  equation  (28)  can  be  rewritten  as 

Aq«  JL  — -  (  ^  \  j  (29) 

c  at  U*/ 

or  using  the  first  order  effect  q  -  KyT  in  equation  (2^) 

Aq  «  T  -fr1-  '  (30) 

n  t 

We  thus  have  a  physical  understanding  of  the  origin  of  the  relaxation  term  in 
equation  (25).  Weymann^1  arrives  at  the  relaxation  time  modified  heat 
conduction  equation  by  considering  the  linear  random  walk  problem. 

To  now  we  have  been  only  considering  what  may  be  called  phonon 
second  sound.  It  may  be  possible  that  electrons  and  holes  in  metals32  and 
semiconductors'^  will  support  similar  macroscopic  thermal  oscillations. 


Let  us  now  consider  equation  (27)  in  some  detail.  Whether  or  not 
finite  velocity  propagating  solutions  are  allowed  depends  upon  the  relative 
magnitudes  of  the  two  terms  having  derivatives  with  respect  to  the  time. 
Clearly  if  the  second  derivative  term  dominates  the  first  derivative  term 
equation  (27)  can  be  approximated  by  the  wave  equation.  Ji  T  is  of  the  form 
T  =  Tq  (r)  exp  (icot),  then  wavelike  solutions  hold  in  the  limits  a)  and/or 
r— ■»,  or  collectively  urr  ^  1* 

To  further  estimate  the  values  of  co  and  r  necessary  for  wavelike 
solutions  let  us  write  equation  (27)  in  the  form 

\  LX  +  J_  XL  -  V?'T  =  0,  (31) 

3 12  x  a  t 

where  v  =  (  — -  /  is  a  velocity,  and  x  -  —  is  the  thermal  diffusivity 

\  /  pC 

H  m 

with  Cm  the  specific  heat  (heat  capacity  per  unit  mass).  The  condition  for 
waveiike  solutions  then  becomes 

vL_  <  1.  (32) 

cox 


Let  us  apply  equation  (32)  to  an  ordinary  material  such  as  room 
temperature  aluminum.  For  v  we  take  the  velocity  of  ordinary  (i.  e.  ,  first) 
sound,  say  6  x  105  cm/sec.  and  X  =  1  cmVsec.  Thus  for  room  temperature 
aluminum  we  expect  that  w>  3.  6  x  1011  radians/sec,  or  a  frequency  greater 
than  5  x  101®  cycles/sec,  is  necessary  for  thermal  wave  propagation.  The 
relaxation  timer  is  actually  the  time  for  energy  or  momentum  loss  (i.  e. ,  the 
relaxation  time  for  the  so  called  "Umklapp"  processes^).  Since  we  must 
have  local  thermal  equilibrium,  we  must  also  satisfy  the  requirement  cor^  <  1 
In  order  to  propagate  a  thermal  wave.  Here  is  the  relaxation  time  for 

momentum  conserving  phonon  collisions  (i.  e. ,  the  so  called  "normal"  processes*^4). 
Thus,  there  is  a  window^'  for  thermal  wave  propagation  gi/en  by 

1  <  oo  ^  1  mm 


1H 


An  appliejd  thermal  frequency  of  approximately  6  x  1011  radians  per  second 
falls  within  the  window  given  by  equation  (33). 

Because  of  the  above  results  we  are  led  to  believe  that  a  high 
energy  laser  pulse,  impinging  on  a  slab  of  aerospace  material,  will  result 
in  energy  leaving  the  skin  depth  material  via  a  propagating  thermal  pulse. 
Such  a  thermal  pulse  would  be  in  addition  to  energy  also  leaving  via  the 
propagating  mass  discontinuity  which  results  from  the  pressure  gradient 
accompanying  the  energy  deposition  in  the  skin  depth  area.  The  energy 
carried  o:‘f  by  the  thermal  pulse  relative  to  that  of  the  propagating  mass 
discontini  ity  is  a  question  yet  to  be  determined,  but  it  is  quite  obvious  that 
the  thermal  pulse  is  not  to  be  neglected  off  hand. 


I 
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ID.  MASS  DENSITY  SHOCKS 


Propagating  mass  density  discontinuities  are  generally  what  people 
have  in  mind  when  they  speak  of  "shock  waves".  The  general  theory  of 
mass  shocks  has  long  since  been  written  down  (e.  g. ,  see  references  7,  21, 
and  36)  and  understood.  As  distinct  from  the  general  theory,  the  state-of-the- 
art  of  mass  shocks  is  concerned  with  effects  which  occur  xn  special  classes 
oi  materials  and  under  special  conditions.  The  question  with  respect  to  these 
effects  is  to  arrive  at  a  model  which  will  explain  the  experimental  observations. 

In  terms  of  state-of-the-art,  this  report  treats  mass  shocks  in  porous 
solids,  various  mass  shock  electromechanical  effects,  and  mass  shocks  in 
energetic  materials  (i.  e.  ,  explosives).  As  areas  of  investigation  they  represent 
a  large  part  of  the  current,  interest  in  "shock"  wave  physics.  For  the  re¬ 
mainder  of  this  report,  by  shock  we  shall  mean  a  mass  shock. 

There  are  three  ways  to  cause  a  shock  wave  to  occur  in  a  material. 

An  already  existing  shock  in  one  medium  can  be  propagated  into  a  second 
material  simply  by  having  the  two  materials  be  in  contact.  Solids  can  be 
caused  to  impact  each  other  at  high  velocity  resulting  in  shock  waves  propagating 
out  from  the  point  of  impact  Lastly,  Ihermal  energy  can  be  stored  in  a  material 
in  a  spatially  inhomogeneous  manner,  since  thermal  energy  is  equivalent  to  a 
pressure,  and  since  dJP  is  equivalent  to  a  force  per  unit  volume  the  inhomogeneous 
dx 

energy  deposition  combined  with  >  nonlinear  equation  of  state  will  result  in  a 
propagating  shock. 

Porous  Solids  -  Any  Materiel  which  has  a  mass  density  less  than  the  maximum 
possible  equilibrium  mass  density  at  a  given  pressure  and  temperature  is  called 
a  porous  solid*.  Examples  er  nuch  materials  are  foams  (e.  g. ,  polyurethane), 
less  than  normal  density  explosives  and  porou.-,  metals. 

Porous  materials  are  of  both  engineering  and  academic  interest. 
Polyurethane  foams  are  utilized  as  shock  mitigators  by  having  the  shock  energy 
absorbed  in  the  process  ot  compacting  the  foam  to  its  normal  solid  density. 

Porous  metals  have  engineering  applications  but  are  also  of  extreme  academic 
importance.  If  a  solid  of  normal  density  is  shock-loaded,  its  final  state  lies 


in  a  single  curve  (the  Hugoniot  curve)  in  the  PV  plane.  On  the  other  hand, 
when  a  solid  such  as  pv  ,vuo  tungsten  is  shock-loaded,  the  final  state 
reached  will  reside  on  one  of  a  continuum  of  Hugoniot  curves,  depending 
upon  the  initial  porosity;  the  irreversible  processes  which  occur  during 
compacting  yield  PVT  states  for  the  compacted  solid  different,  from  those 
occurring  in  the  shock-loaded  material  having  initial  normal  density.  Thus, 
by  shock-loading  an  initially  porous  material  to  beyond  complete  compaction, 
it  is  possible  to  investigate  bread  regions  in  PVT  space.  This  academic 
property  of  porous  materials  was  first  pointed  out  by  Zeldovich37. 

Another  important  application  oi  porous  materials  is  in  the  area 
of  explosives.  It  is  well  known  that  the  initiation  sensitivity  of  explosives  to 
mechanical  shock  depends  upon  the  final  deusity  to  which  a  granular  explosive 
may  be  pressed^.  This  area  also  has  both  academic  and  engineering 
applications. 

Soviet  and  Western  interest  in  inert  porous  materials,  judging  from 
the  open  literature,  has  been  in  different  pressure  regimes.  Western  interest 
has  been  strong  in  the  area  of  shock  propagation  in  the  pressure  range  at  which 
compaction  occurs.  On  the  other  hand,  the  Russians  appear  to  have  their 
main  effort  in  the  pressure  range  beycad  which  compaction  has  already  occurred. 
Among  the  interesting  ingredients  of  high-pressure  physics  which  the  Russians 
have  studied  via  porous  materials  is  the  electronic  contribution  to  the  risen 
constant36’  33. 


Presently  the  major  unkiowr.  in  the  area  of  porous  solid  physics 
concerns  the  mechanism  of  the  compaction  process.  In  a  typical  experiment 
a  one-dimensional  strain  shock  (no  net  local  particle  displacement  per¬ 
pendicular  to  the  direction  of  shock  propagation)  the  free  surface  velocity  and 
an  average  shock  velocity  are  measured.  Then,  by  using  Rankine- Hugoniot  type 
conservation  equations7’  46  (for  mass,  momentum,  and  energy  flows)  the 
pressure  and  energy  density  may  be  determined'*6*  in  the  shocked  regime. 
Alternatively,  the  porous  material  may  be  backed  up  with  a  quartz  (or  other 
piezoelectric)  transducer,  to  yield  pressure  and  panicle  velocity Again, 
the  conservation  equations  are  used  to  yield  shock  velocity  and  energy  density. 


Thvs,  experimentally  it  is  possible  to  catagorize  a  porous  solid 
according  to  its  gross  macroscopic  behavior  in  a  shock  environment.  The 
motivation  for  fhe  axuenments  usually  is  to  arrive  at  or  to  verify  a  pro¬ 
posed  equal  .on  of  eta  to  in  the  dynamic  regime.  Such  an  equation  of  state  is 
generally  of  the  form 


P  -  P  * 


(34) 


where p  is  the  actual  locai  moss  density,  p  is  the  mass  dens' ^  “  v  the 
material  if  it  voidless  but  o*.  erwise  iu  the  same  studo  as  the  j  solid 
under  consideration,  and  £  is  the  internal  energy  density.  me  depen¬ 

dence  is  meant  to  include  strain  "•-ite  effects. 


The  type  of  porous  material  equation  of  ...tate  modeling  which  ha3  so 
far  occurred  centers  on  choosing  particular  functional  forms  tor  conceptual 
physical  models)  which  correspond  to  equities  <34).  One  can  think  of  the 
existing  modeling  as  belonging  to  three  different  classes. 

a.  The  locking  solid  model  !3.  In  thi.i  model  the  P  versus  p 
depv  ’snee  is  basically  taken  as  shown  in  figure  4.  Variations  m  tlfe  model  are 
achieved  by  allowing  the  compacted  material  to  be  more  or  less  compressible, 
and  by  allowing  the  yield  strength  tc  be  more  or  less  significant. 
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Figure  4  -  The  lockii.r  mod**!  tor  porous  solid  compression.  Y  is  the 

yield,  sties-  and  Vs  is  the  volume  of  the  original  undeformed 
rnat*'**iai  we  it  .otdlesj 
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b.  The  P-«  model  ,  Whereas  the  locking  modal,  as  shown 
in  figure  4,  has  «  =  P  s  independent  of  pressure  in  the  compaction  i  egime, 

P 

the  P-«  model  gives  «  an  explicit  pressure  dependence.  In  particular.,  if  one 
assumes  a  polynominal  expansion  for  «  ’  ' 

oc  =  «  +  <c  P  +  o:  +  a  p^  +  .  .  .  ,  (35) 

o  1  !2  3 

then  suitable  boundary  conditions  on  «  and  d«  (evaluated  at  the  yield  and  fully 

dP 

compacted  points  in  the  P-p  plane)  are  ufficient44  to  determine  ,  Wj  ,  «  , 
and  «g  . 

1  \  c.  The  plate  -  air  gap  model45.  In  this  model  the  porous 

material  is'  viewed  as  plates  of  solid  material  (the  solid  matrix  of  the  porous 
medium)  separated  by  air  gaps  (corresponding  to  the  voids).  A  1-D  strain 
shock  propagates  via  the  plates  being  transported  across  the  air  gaps,  in 
effect  collapsing  the  voids,  anu  thus  resulting  in  a  voidless  material  behind  the 
shock  front 

The  locking  solid  model  and  the  P-«  model  have  in  common  that  they 
are  .  asicatly  empirical  equation  of  state  models.  They  are  not  based  upon  the 
microscopic  details  of  the  porous  solid  structure  in  the  spirit  of  arriving  at 
an  equation  of  state  of  a  crystalline  solid  via  a  knowledge  of  interatomic 
potential®®.  On  the  other  hand  the  plate-air  gap  model,  including  more 
sophisticated  versions  of  it4®,  is  based  upon  a  macroscopic  equivalent  of  a 
microscopic  model. 

While  each  of  the  models  mentioned  above  is  capable  of  predicting 
shock  velocity  and  attenuation,  none  of  them  are  capable  of  predicting  either 
the  rise  time  (  and  thus  shape)  of  an  elastic  precursor  or  the  main  compaction 
wave.  While  an  attempt  has  been  made47  to  include  a  stress  relaxation  term 
within  the  P-«  model,  to  now  the  approach  has  only  been  empirical  and  the 
fittod  relaxation  times  have  no  obvious  relation  with  z'espect  to  the  actual 
details  of  void  collapse4®. 
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In  order  to  really  understand  the  compaction  process  in  porous  c 
materials  it  seems  necessary  to  first  consider  the  detailed  deformation 
mechanics  of  individual  voids  in  a  representative  matnx  material,  and 
then  to  develop  suitable  methods  of  averaging  over  the  entire  material. 

With  such  an  approach  the  necessary  time  dependent  constitutive  equations 
(i.  e. ,  explicit  strain  rate  dependence),  the  shock  velocity,  attenuation, 
and  wave  shapes  would  automatically  result.  The  approach  has  already  been 
successful  for  shock  calculations  in  solids  having  dislocations'^. 

Energy  deposition  experiments  (one  of  the  three  methods  of  shock 
generation  previously  mentioned)  in  porous  solids  have  shown  some  interesting 
effects.  It  has  been  found  experimentally^1  ^  that  the  Gruneisen  constant 
of  porous  solids  is  very  much  dependent  upon  the  degree  of  porousity.  The 
Gruneisen  constant,  F,  is  defined  thermodynamically  as 


p  tae  A 


where  the  derivative  is  to  be  evaluated  at  constant  density.  Table  II  below 
gives  some  of  the  experimental  data  for  porous  PETN  obtained  by  using  an 
electron  beam  to  deposit  energy  in  the  material.  Although  some  of  the  data 

shown  in  Table  II  have  very  large  errors,  the  trend  as  a  function  of  p  /a 

o  30 

clearly  discernible. 

TABLE  IIa 

ENERGY  DEPOSITION  PROPERTIES  OF  PETN 


/  b 
p  //B 

Density 

Sound  Speeds 

Effective 

o  80 

Era/crn^ 

cm/sec  x  10~a 

r 

0.  95 

1.67 

2.  8 

1.2 

0.  90 

1.59 

2. 4  +  0. 1 

0.51 

0.  87 

1.54 

1.  8  +  0.  3 

0.  15 

0.  84 

1.  48 

1.7  +0.3 

(0,  07  to  0.  23) 

aThis  table  was  taken,  with  minor  changes,  from  reference  42. 

bp  is  the  density  of  the  undeformed  porous  PETN,  while  p  is  the  density 
o  so 

of  undeformed  voidless  PETN. 
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The  experiments  of  Shea  et  al*  resulting  in  Table  II  are  interesting 
in  that  the  electron  beam  was  depositing  energy  for  a  longer  period  of  time 
than  it  takes  a  sound  wave  to  propagate  across  a  typical  particle  of  the  porous 
PETN  material.  Consequently,  the  T  measured  was  not  the  true  constant 
volume  Gruneisen  parameter  as  defined  by  equation  (36).  Exactly  what 
detailed  microscopic  physics  corresponds  to  the  measured  r  is  a  very 
difficult^®  question  in  porous  solid  mechanics  which  remains  unanswered. 
While  the  measured  "effective  F"  is  a  worthwhile  experimental  variable  in 
that  it  characterizes  the  gross  behavior  of  the  material,  unforturnately  it 
does  cot  allow  for  distinguishing  between  all  possible  porous  materials  as  a 
function  of  material  parameters. 

Mazella  et  al5®  have  proposed  a  Gruneisen  parameter  for  a  porous 
material  given  by 

c2 

r=  ~T  rs  ,  (37) 

c 

s 

where  Fg  is  the  parameter  for  the  voidless  undeformed  solid,  with  cs  being 

the  propagation  velocity  corresponding  to  a  pure  bulk  modulus  wave51ain  the 
voidless  undeformed  material,  c  is  the  corresponding  velocity  for  the  porous 
solid.  Equation  (37)  follows  from  the  assumed50  form  for  the  equation  of 
state  in  differential  form 


(38) 


where  K  denotes  bulk  modules.  Equation  (38)  differs  from  the  equation  of 
state  of  a  voidless  material  by  the  presence  of  the  factor _K_  .  Physically 

K 

s 

the  factor  K_  spreads  the  energy  which  the  electron  beam  has  deposited  in 
K 

a 

the  solid  particles  of  the  porous  material  over  the  entire  solid.  To  see  this 
consider 


K  -  V 


SP 

av 


ap 

av 


(ap) 


(39) 
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where  (AP)*  indicates  an  equal  differential  pressure  to  porous  and  voidless 
media  having  the  same  solid  composition.  By  combining  the  energy  density 
term  of  equation  (38), 


AP  = 


<AP)g  . 


(40) 


with  equation  (39)  we  find 

<&)  =  *-£-  (API  or  e  (AP)  =  c  (AP)  ,  (41) 

v  8  bo 

where  €g  is  the  strain  in  the  voidless  material  for  the  increment}?,!  (AP)*. 

The  second  equation  (41)  is  the  energy  statement  alluded  to  above. 

It  thus  appears  that  equation  (37)  is  valid  for  those  experimental 
situations  where  the  energy  deposition  is  so  long  in  time  that  the  porous 
medium  is  able  to  continuously  adjust  its  density  on  a  microscopic  scale 
without  effecting  the  macroscopic  density.  In  some  experimental  cases 
equation  (37)  is  found  to  hold,  and  in  others  not  to  hold50.  Unfortunately  we 
are  not  given  sufficient  information  about  the  materials  examined  (e.  g. , 
particle  size)  to  judge  the  theoretical  applicability  of  equation  (37)  via 
equation  (41).  Even  though  we  expect  a  time  dependent  compaction  process, 
equation  (37)  does  not  contain  a  relaxation  time. 

Electromechanical  Efteots  -  Experimentally  it  is  well  known  that  electrical 
signals  originate  in  many  materials  when  a  shock  wave  is  caused  to  propagate 
through  the  material.  No  external  source  of  electric  power  is  necessary  for 
the  existence  of  such  signals,  a  typical  experimental  arrangement  for  the 
effect  in  question  is  shown  in  figure  5. 

An  obvious  example  of  a  specimen  material  is  quartz,  because  of  its 
piezoelectric  character.  The  effect  in  quartz  is  so  *vell  understood  that  quartz 
is  used  as  a  ga^ge^^to  measure  shock  amplitudes  in  many  experiments.  On 
the  other  hand,  many  materials  which  are  not  piezoelectric  (some  have  highly 
symmetric  crystal  structures)  also  yield  el  ictrical  signals  upon  shock  loading. 
The  alkali-halides5^,  and  germanium  and  silicon5^  are  examples  of  non- 
pi  ozoelectrics  in  which  the  effect  is  seen.  Indeed,  even  plastics54  yield 
signals. 


A  =  Ammeter 


SZ&ifVFmir! 


V  =  Voltmeter 
R  =  Load  Resistor 

Ij 


FIGURE  5  -  Typical  Experimental  Arrangement  for  Measuring 
Shock  Induced  Electrical  Effects 
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There  also  is  an  interesting  related  effect  known55  as  the  "anomalous 
thermoelectric  effect".  Here  a  thermoelectric  junction  is  shock  loaaed  with 
the  shock  propagating  perpendicular  to  the  plane  of  the  junction.  The  observed 
electrical  signal  is  approximately  an  order  of  magnitude  larger  than  that 
associated  with  the  thermoelectric  voltage  which  would  be  associated  with  the 
compression  caused  adiabatic  temperature  rise  in  the  junction  materials. 

To  date  explanations  for  the  effect  in  plastics  do  not  exist.  A  crude 
qualitative  theory  involving  dislocation  motion  has  been  put  forward55  to 
explain  the  shock  induced  voltages  in  the  centrosymmetric  alkali-halides,  and 
first  order  transport  theory  has  been  applied  to  germanium5?,  silicon,  and 
metallic58  band  structures  as  first  attempts  which  show  the  theoretical 
existence  of  the  effects  ii  those  materials.  The  theory  for  germanium  and  silicon 
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is  essentially  an  extension  of  the  so-called  "acousto-electric"  effect®® 
wherein  experimentally  the  effect  is  seen  via  acoustic  stimulation.  The 
semiconductor  predictions  are  in  order  of  magnitude  agreement  with  the 
experiments  of  Mineev  et  al®3,  while  the  theoretical  results  for  the  metallic 
case  are  still  open  to  question.  For  the  metallic  case  Harrises  predicts 
approximately  10  microvolts  per  100  kilobars  of  impressed  shock,  a 
prediction  orders  of  magnitude  smaller  than  that  observed  with  shocked 
metallic  thermocouple  junctions®®. 

For  the  metallic  thermocouple  junction,  Migault  and  Jacquesson®® 
have  attempted  to  explain  the  anomal  is  thermoelectric  power  by  invoking 
a  pressure  induced  force  on  the  conduction  electrons.  Their  treatment, 
however,  does  not  take  into  account  relaxation  time  effects  within  and/or 
behind  the  shock  front  -  a  physics  point  which  was  found  to  be  of  paramount 
importance  (i.  e. ,  zero  voltage  if  it  was  neglected)  for  the  case  of  a  single 
metal®®.  The  Migault  and  Jacquesson  results  are  in  order  of  magnitude 
agreement  with  experiment. 

Directed  to  the  anomalous  thermoelectric  power  problem,  Conze 
et  al®1*  ®2»  ®3  have  done  some  very  general  work  concerning  electrons  and 
electron-phonon  interactions  in  a  shocked  solid.  Directed  towards  the  "shock 
polarization"  of  the  alkali-halides,  Harris®4  has  considered  the  possibility  of 
the  elastic  orientation  of  the  Hydroxide-ion  impurities  which  are  generally 
found  as  unwanted  impurities  in  the  as-grown  crystals.  Such  a  theory  would 
require  a  nonsymmetric  stress  tensor,  and  make  for  a  very  interesting 
theoretical  mechanics  problem  (in  spite  of  experimental  evidence  indicating  that 
the  effect  is  unimportant®®).  In  addition  to  the  work  of  Harris®7  for  the 
semiconductor  problem,  Horie66  has  published  work  based  upon  a  phonon- 
drag  type  of  effect  which  is  in  the  right  direction  theoretically.  Horie  believes 
that  his  results  indicate  no  measureable  signals  at  room  temperature,  in 
contradiction  to  experiments®3  where  results  are  indeed  seen. 

The  question  of  shock  induced  electrical  effects  ir.  solids  is  wide  open. 

As  of  this  writing  there  are  not  any  microscopic  theories  capable  of  predicting 
electrical  voltages  for  arbitrarily  large  shock  amplitude.  The  best  which  can 
be  said  about  the  present  situation  is  that  some  first  or  second  order  theories 
allow  one  to  accept  the  existence  of  the  shock  induced  voltages  on  a  theoretical 
basis.  This  area  of  shock  physics  requires  a  great  deal  more  detailed  theoretical 
and  experimental  work  to  be  done  before  one  could  feel  remotely  happy  about 
the  situation. 


Energetic  Materials  -  Porous  materials  and  electromechanical  effects  were 
considered  first  in  this  chapter  because  both  subject  areas  are  considered  to 
be  important  to  the  physics  of  initiation  of  energetic  materials  (explosives). 

The  fundamental  microscopic  processes  which  describe  initiation 
and  detonation,  in  terms  of  the  parameters  which  define  an  undetonated  chemical 
explosive  material  are  not  understood.  Indeed  the  thermonuclear  detonation 
wave  structure®7  is  much  hotter  understood  than  that  of  the  chemical  detonation 
wave  structure.  The  reason  for  the  above  surprising  statement  is  simple; 
nuclear  scattering  cross  sections  and  processes  (millions  of  electron  volts)  are 
better  understood  then  chemical  "scattering"  oroea  sections  and  processes 
(electron  volts).  Furthermore,  in  the  electron  volt  regime  the  energies  con¬ 
tributed  by  solid  state  coherent  effects  are  expected  to  be  of  importance  while  no 
such  statement  can  be  made  for  the  Mev  regime.  In  other  words  extended 
intramolecular  as  well  as  intermolecular  effects  are  expected  to  be  important 
in  the  processes  which  characterize  solid  chemical  explosives. 

Semi-recent  literature  has  contained  experimental®®  and  theoretical®®*  7® 
hints  of  electromechanical  phenomena  playing  a  role  in  the  detonation  process. 
Unfortunately  the  type  of  experiments  so  far  undertaken  are  not  the  rigorous 
solid  state  experiments  that  a  theoretical  physicist  would  like  to  work  from. 

Rather  the  experiments  point  only  to  the  existence  of  possible  ionization  phenomena 
associated  with  initiation  to  detonation.  The  theoretical  work  has  to  now  only 
dealt  with  individual  mechanisms  (e.  g. ,  graded  band  gap  effects7®)  without 
completely  carrying  the  calculations  through  to  detonation  physics.  It  is 
interesting,  however,  that  Williams7®  in  effect  suggests  that  an  electric  dipole 
layer  may  be  associated  with  the  detonation  front  while  Mineev  et  al®®  and  Harris®7 
predict  a  dipole  layer  to  be  associated  with  the  shock  front  in  an  inert  material. 

Other  approaches  to  the  question  of  initiation  are  also  taking  place. 

For  example,  various  groups  are  studying  the  crystal  struct  ,re  of  solid 
explosives  via  neutron  diffraction  techniques’'1,  while  some  workers  are 
considering  the  detailed  micromechanics  of  bond  breakage  under  a  plate  impact 
type  of  environment72.  Aside  from  giving  detailed  crystallographic  information, 
there  is  the  possibility  that  neutron  diffraction  experiments  will  lead  to  the 
discovery  of  a  "soft  mode"  and  a  consequent  lattice  instability73.  Of  all  the  non- 
tbermal  approaches  to  the  question  of  initiation  to  detonation,  it  would  appear 
that  this  crystallographic  approach  is  the  most  rigorous  and  complete  so  far 
undertaken . 
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To  date  the  most  successful  attack  on  the  initiation  problem  appears 
to  bf  via  the  macroscopic  thermal  route,  sometimes  called  "hot  spot"  theory. 
Trot'  et  al*1  have  shown  that  an  Ahrenius  term  such  as  exp(-A/RT),  when 
combined  with  adiabatic  heating  effects  resulting  from  plate  impact  caused 
compression,  is  sufficient  to  explain  the  observed  nonelectrical  phenomena 
associated  with  homogeneous  and  heterogeneous  explosives.  The  success  goes 
so  far  as  to  enable  the  calculation  of  a  reaction  or  incubation  time  for  detonation 
to  occur  at  the  impacted  surface  in  a  homogeneous  explosive. 

The  question  of  the  microphysics  pertinent  to  detonation  phenomena,  in 
terms  of  the  parameters  which  are  important  to  solid  state  physics,  is  wide  open. 


VL  SUMMARY 


We  have  surveyed  the  field  of  shock  wave  physics  in  solids.  In 
so  doing  we  have  generalized  the  concept  of  a  shock  wave  to  include  thermal, 
electromagnetic,  as  well  as  the  more  generally  understood  mass  density 
shocks. 


We  have  briefly  examined  the  state-of-the-art  with  respect  to 
understanding  the  compaction  behaviour  of  porous  materials,  the  micro¬ 
scopic  explanation  of  electromechanical  effects  (shock  polarization)  in  a 
variety  of  media,  and  the  state  of  our  understanding  cf  the  physics  of  the 
detonation  process.  Porous  materials,  electromechanical  phenomena,  and 
detonation  physics  represent  wide  open  areas  for  future  investigation. 
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